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We show that weak off-diagonal disorder in degenerate ground state conjugated polymers results 
in a finite density of randomly positioned kinks (solitons and antisolitons) in the lattice dimerization. 
For realistic values of the disorder, these kinks should clearly show up in the optical and magnetic 
properties. 

PACS numbers: 71.20.Rv, 63.20.Kr, 71.23.-k 

For over thirty years, conjugated polymers have provided a rich field of fundamental and technological research 
jjj. They belong to a large class of quasi-one-dimensional Peierls materials, in which the lattice distorts due to its 
interaction with itinerant electrons Much attention has been devoted to trans-polyacetylene, which has a doubly 
degenerate ground state. The degeneracy allows for topological excitations, solitons, which are kinks in the lattice 
dimerization accompanied by a local distortion of the electron density H. Solitons provide a simple explanation of 
the unusual properties of doped trans-polyacetylene |]||. They also feature in quantum lattice fluctuations, which 
lead to subgap optical absorption Q and may be responsible for an enhanced nonlinear optical response |^| . Solitons 
were first found in the free-electron SSH model |5]||. Their topological nature, however, allows them to survive also 
in the presence of on-site Coulomb repulsion (Peierls-Hubbard model) . 

Although real polymer materials contain many types of defects, most theoretical studies have been restricted to 
perfect chains. Due to the quasi-one-dimensional nature of polymers, however, any disorder will play a major role in 
determining their properties. In this Letter, we study the effects of weak disorder in the electron hopping amplitudes in 
polyacetylene chains. Such (off-diagonal) disorder arises from random chain twists that diminish the overlap between 
the 7r-orbitals of neighboring carbon atoms. In previous analytical studies, off-diagonal disorder was modeled by small 
fluctuations of the Peierls order parameter around an average value |^||. Numerical calculations showed, however, 
that if the lattice configuration of a weakly disordered chain described by the SSH model is allowed to relax, it may 
contain solitons, i.e., large (nonperturbative) variations of the order parameter JTo| , pT[ . This implies that the lattice 
response to the disorder cannot be neglected. An analytical theory for this phenomenon and an estimate (either 
numerical or analytical) for the density of thus created kinks, have sofar not been given. In this Letter, we address 
these issues and we show that, in fact, arbitrarily small off-diagonal disorder leads to a finite density of solitons. Our 
arguments also hold in the presence of electron correlations. 

We consider a half-filled chain, in which itinerant (it) electrons interact with the classical carbon lattice. We assume 
that the chain Hamiltonian contains three contributions: (i) the electron hopping between adjacent carbon atoms, 
(ii) the elastic energy of the lattice, and (iii) the Coulomb repulsion between the electrons. The hopping amplitude 
has the form, 

t m ,m+l = to + ce{u m — u m+ i )+St 

rn.rn-\-l ■ (i) 

Here the first term is the bare amplitude and the second term describes the electron-phonon interaction, where a is 
the coupling constant and u m denotes the displacement of the m'th carbon atom from its uniform-lattice position. 
Finally, the third term is a random contribution, which we include to describe the conformational disorder. We assume 
that the fluctuations <5t m ,m+i are frozen ("quenched" disorder) and that they are independent for different links. We 
will restrict ourselves to zero-temperature properties. 

As we shall shortly see, for weak disorder the spatial extension of a disorder fluctuation necessary to create a kink 
is much larger than one lattice unit, which allows us to use a continuum model. Then Eq.(|l]) is replaced by, 

A(x) = A lat (x) + r)(x) , (2) 

where the lattice dimerization A/ at (2met) = a(u2m-i — 2u2m + U2m+i) (a is the lattice constant) describes the 
alternating part of the hopping amplitude determined by the shifts u m of the carbon atoms (cf. the continuum 



version of the SSH model |12|) and ry(2ma) = <5i2m-i,2m — ^2m,2m+i describes the disorder part with a Gaussian 
correlator, 

( V (xHy)) = A6(x - y) . (3) 
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Note, that while chain twists always decrease the hopping amplitudes (5t m . m+ i < 0), T)(x) can be both positive and 
negative, as it is the alternating part of the fluctuations. 

It should now be realized that in a consistent treatment of disorder, A; Qt (x) implicitly depends on the rj(x). The 
reason is that at zero temperature Ai at (x) for a particular disorder realization is found by minimizing the total 
(electronic plus lattice) energy. The total energy reads, 



E = E el [A(x)j + — — / dx Af at (x) , (4) 
ttXv f J 

where A = ^ K is the dimensionless electron-phonon coupling constant (K is the spring constant), vp = 2ato is the 
bare value of Fermi velocity, and we set h = 1. The electron energy depends on A(x) defined by Eq.(||), which includes 
the changes of the hopping amplitudes due to both the lattice distortion and the disorder. We stress that in what 
follows we will not need an explicit expression for the electron energy, which enables us to include the electron-electron 
interaction. We will assume, however, that in the absence of disorder the total energy of a half-filled chain has two 
minima, both corresponding to a uniform dimerization, A(x) = ±Aq. We will also assume that the model admits 
soliton solutions, i.e., apart from the above minima, the total energy has an infinite number of extrema A^(x\z) 
corresponding to the presence of N kinks (solitons and antisolitons interpolating between — Ao and +Ao and vice 
versa), whose positions are described by the iV-dimensional vector z = (zi, z%, . . . , z^). Both the SSH and the Peierls- 
Hubbard model have these properties. If the average separation between neighboring kinks is large compared to the 
correlation length £o = vpj Ao, the configuration with N kinks has energy, 

E N =E + NE S , (5) 

where E s is the energy needed to create a single soliton. The kinks can be either charged and spinless or neutral with 
spin ~ |J. In the SSH model, the two types have the same energy, while in the presence of Coulomb repulsion the 
neutral soliton is energetically favorable J7],|l3| and in undoped chains will be the only type of kink that occurs at zero 
temperature. 

At weak disorder, \rj(x)\ <g; Ao, linear response theory predicts small fluctuations of the ground state lattice 
configuration around the uniformly dimerized lattice. We will show, however, that for arbitrarily weak disorder the 
ground state configuration cannot be found using linear response theory, as the configuration with small fluctuations 
around a chain containing a number of solitons may have a lower energy than the one with small fluctuations around 
uniform dimerization. 

To first order in ^(x), the correction to the energy of a multikink configuration reads, 



5E N = — dxA N (x\z)r)(x) , (6) 

n\v F J 

where the extremum condition for the ordered configuration Ajy(x\z) was used. Hence, the change of the total energy 
of the uniformly dimerized state (A(x) = Aq) due to disorder reads, 

SE = --^2. fdxr,(x) . (7) 

Consider now, on the other hand, the same chain with an antisoliton at z\ and a soliton at Z2, chosen in such 
a way that the whole disorder fluctuation lies between z\ and Z2- The change of energy due to disorder for this 
configuration equals —5E 0l because between z\ and z 2 A 2 (x\zi, z 2 ) ~ — A . We thus find that in the disordered 
chain, the configuration with a soliton-antisoliton pair is energetically favorable to the uniform configuration if 

dxrj(x) > jXvp . (8) 



Here 7 = which in the case of the SSH model equals 1. It should be stressed that as long as the kink density 
is indeed small compared to l/£,o, the entire effect of electron-electron interactions is contained in the factor 7. We 
note that, no matter how small the disorder is, the fluctuations in Jdxr](x) grow with the chain size, so that for a 
sufficiently long chain the inequality Eq.(||) will certainly be fulfilled. 

The creation of kinks by off-diagonal disorder is illustrated in Fig. 1. The thick line is the ground state order 
parameter A rn = £2m-i,2m ~ ^2m,2m+i f° r a discrete chain of 160 carbon atoms, obtained by numerically minimizing 
the total energy for one particular realization of the disorder. In this example, we neglected the Coulomb repulsion 
and used standard SSH parameters. It is clearly seen that the order parameter fluctuates near a soliton-antisoliton 
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pair configuration. The disorder realization that was used as input is visualized by the thin line, which gives Ao+rj m — 
Ao + St2m-i,2m — ^2m,2m+i (^o ~ 0.7eV). The thin line in fact directly represents the order parameter along the 
chain if we neglect the lattice relaxation. In that case, no lattice kinks are created. 

In order to obtain an estimate for the density of disorder-induced kinks, we consider a disorder fluctuation with 
spatial size I and constant value r\\ . Then, according to Eq. (||) , the threshold for the creation of a soliton-antisoliton 
pair is, 

\m\ = — — ■ (9) 

By requiring the probability density of such a fluctuation, 

1 / , J^„^^2^ ( (i^f) 2 



p = exp ( -— dxri(x) J=cxp[ ) , (10) 



to be of the order of unity, we obtain the average density of kinks: 

1 A 



I (7Ai>i?) 
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(11) 



The condition p ~ 1 ensures that the typical value of $ x+l dx'ri(x') is large enough to stabilize a soliton-antisoliton pair. 
The probability to find such a pair with separation much smaller than the average given by Eq.(|ll]), is suppressed 
(cf. Eq. (|To|) ) as it requires a large disorder fluctuation. On the other hand, lattice configurations in which the average 
distance between the kinks is significantly larger than I given by Eq. ([□]) are statistically less likely. Similar arguments 
were used to estimate the typical domain size for the Ising model in a random magnetic field |I3| . The analogy is not 
accidental: it can be shown that the statistical properties of the disorder-induced kinks are indeed described by the 
latter model. 

We next address the question whether our first-order expansion around multi-kink configurations is a meaningful 
approach. First, we estimate second-order effects. To this end, we again consider the fluctuation with spatial size I 
and magnitude |ry;|. In this case, the ratio S^E/S^E of the second- and first-order corrections to the total energy 
can be expressed in terms of the renormalized frequency of the phonon with wave-vector ir/a. This frequency can be 
calculated in the two limiting cases for the on-site electron-electron interaction U. We then find: 

S^E_tof |(1 -A), for U = (12) 



6<VE I \ \l\ for C7»* 

plus higher order terms in powers of £o/l. Therefore, our approach is valid as long as the average distance between 
the kinks is much larger than the correlation length (dilute gas of kinks) . 

Second, we should stress that, although the energies of the multi-kink configurations are calculated perturbatively, 
the stabilization of kinks by weak disorder is clearly a non-perturbative effect. The reason is that for growing disorder 
the minimal-energy lattice configuration with a finite density of kinks does not evolve continuously from the uniformly 
dimerizcd configuration. In fact, as long as the second-order correction to the energy due to disorder is relatively 
small, each multikink lattice configuration docs not change much in the presence of weak disorder. Nevertheless, 
their energies, even in the linear regime, may change appreciably for sufficiently long chains. This fact is illustrated 
in Fig. 2, where we plot the energies of several multi-kink configurations, with to 4 soliton-antisoliton pairs, as a 
function of the magnitude of the disorder. Here, the shape of the disorder fluctuation is assumed to be fixed, and 
the slopes of the plotted lines depend on the particular choice of the kink positions. From this figure it is clear that 
the increase of the number of soliton-antisoliton pairs in the minimal-energy configuration for a given shape of the 
disorder fluctuation occurs discontinously in the region where the energies are still linear in the disorder. 

Elsewhere, we show in detail that the estimate for the density of disorder-induced kinks Eq.(pT|) is, in fact, the 
exact result in the limit where the average distance between kinks is much larger than the correlation length £o (dilute 
gas of kinks) and much smaller than the length of the chain L. The exact solution was obtained by noticing that 
the disorder averaged free energy of the chain at sufficiently low temperature has the form of a matrix element of 
the Green function describing the relaxation of a spin-i in a random magnetic field. The coordinate along the chain 
plays the role of the (imaginary) time in which the relaxation takes place, while the kinks correspond to spin-flips. 
The Green function was found by solving the corresponding Fokker-Planck equation. 

We now discuss the consequences of the disorder-induced solitons for the optical and magnetic properties. In the 
absence of electron-electron interactions, the finite density of kinks leads to the appearance of a peak in the single- 
particle density of states at e = 0, because each kink carries an electron state with zero energy M. At weak disorder 
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this peak is sharp since the distances between the kinks are large, so that the splitting of the energies of the electron 

states localized near different kinks due to their mutual interaction is very small (~ exp ^— -^j)- The peak in the 

density of states should result in a midgap peak in the optical absorption spectrum. In undoped trans-polyacetylene, 
however, this peak is not observed. This can be ascribed to a small magnitude of disorder, which is, however, unlikely, 
since the average (presumably disorder limited) conjugation length of trans-polyacetylene is known to be of the order 
of several tens of carbon atoms p5[ . 

The on-site Coulomb repulsion decreases the neutral soliton energy (?],[k| (7 < 1), thus increasing the kink density 
(Eq.(fflJ)). However, at the same time, the Coulomb interaction shifts the absorption peak resulting from the neutral 
kinks towards higher energies, where it may merge with the peak resulting from interband transitions. This may 
explain the absence of a clear midgap peak in the absorption spectrum, but it does not reconcile the theory with 
magnetic measurements. The neutral disorder-induced kinks have spin 1/2 and should contribute to the magnetic 
susceptibility. However, the observed Curie susceptibility for undoped trans-polyacetylene corresponds to only one 
free spin per 3000 carbon atoms Jl6| , which is much smaller than the average density of kinks given by Eq. ([ll|) for 
the disorder magnitude used in [EI to fit the absorption spectrum. 

This contradiction between theory and experiment clearly poses questions as to the applicability of the SSH or 
the Peierls-Hubbard models to trans-polyacetylene. Here we note several factors that may solve this problem. First, 
interchain interactions counteract the formation of free kinks, binding them into pairs. Second, quantum lattice 
motion induces an effective interaction between the spins of neighbouring solitons, resulting in a spin-liquid with a 
singlet ground state. Finally, the free spins of solitons can be bound by chemical impurities. Further research is 
needed to see whether the quantitative effects of these factors are strong enough to account for the strong reduction 
in the number of observed free spins. 

In summary, we have shown that in Peierls chains arbitrarily weak off-diagonal disorder induces a finite density of 
neutral solitons. As we discussed, these solitons strongly affect the observable properties of the chains. We believe, 
that they also play an important role in the phase transition of disordered materials into the Peierls state. 
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FIGURE CAPTIONS 

FIG. 1. Numerically obtained order parameter along a polyacetylene chain of 80 pairs of carbon atoms for one 
particular (but typical) realization of off-diagonal disorder. The thin line does not account for lattice relaxation, the 
thick line does. 

FIG. 2. Illustration of the energies of several multi-kink configurations as a function of the magnitude of the disorder 
fluctuation. From bottom to top (at zero disorder), the curves correspond to chains with to 4 soliton-antisoliton 
pairs. 
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